Phase structure of a quantized chiral soliton on S"^ 
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A quantization of a breathing motion of a rotating chiral soliton on 5''^ is performed 
^ I in terms of a family of trial functions for a profile function of the hegdehog ansatz. 
^ ■ We determine eigenenergies of the quantized skyrmion by solving the Schrodinger 

m 



equation of the breathing mode for several lower spin and isospin states varying the 



i-C I Skyrme term constants e. When radius is smaller than Ijej^^, where j^^ is the pion 
P-i! 

Q_|' decay constant, we always obtain a conformal map solution as the lowest eigenenergy 
(D ! 

^ . state. In the conformal map case, allowed states have only symmetric or anti-symmetric 



wave function under inversion of a dynamical variable describing the breathing mode. As 
the radius increases the energy splitting between the symmetric and anti-symmetric 
states rapidly decreases and two states become completely degenerate state. When the 
radius larger than S/c/tt, for the small Skyrme term constant e the lowest eigenenergy 
states are obtained with the profile function given by an arccosine form which is almost 
the same to those of usual Ft? skyrmion. When the effects of the Skyrme term are weak, 
i.e. large e, the lowest energy states are obtained by the profile function of conformal 
map, which correspond to the "frozen states" for the F? skyrmion as the limit of 
radius oo. 
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§1. Introduction 

A minimal extension of the nonlinear chiral Lagrangian can admit soliton 
solutions called skyrmions which have been considered as baryons, and the Skyrme 
model has been shown to provide a reasonable phenomenology for static properties of 
baryons.^-* Since Klebanov applied the skyrmion to the extended system of nuclear mat- 
ter with a finite baryon density, a variety of numerical studies of the periodic arrays of 
the cubic lattice skyrmion in the flat space have been made.^)~^^ These approaches, 
however, are technically so complicated and require tedious numerical works due to the 
huge number of dense mesh points. The other approach is initiated by Manton where 
the flat space is replaced by the S-sphere.^^'^^ In this approach a single baryon on 
a 3-sphere S^{p) with radius p is just an approximation of baryonic matter with an 
average baryon density of l/27r^p'^ and provides very good approximation of the generic 
features including the phase transition of the R^ lattice skyrmion approaches. 

The most striking fact that is found from these studies is existence of two phases. 
In the skyrmion case, one is a phase where both baryon number and energy densities 
are well localized around the north or south pole of corresponding to the low baryon 
density matter. The other is a phase where these densities are much more homoge- 
neously distributed and corresponds to the high baryon density matter. It is expected 
that transitions between the localized states and the homogeneous states will be closely 
related to the breathing motion of skyrmion. So far, as we know, no studies of the 
skyrmion on quantum level including breathing mode have been made. Recently, it has 
been shown by Yang and one of the author^^ that for a certain region of the Skyrme 
term constant the generic features of the R^ skyrmion on quantum level is drastically 
changed from those on classical level and that the rotational and breathing modes are 
completely frozen. Therefore, it is very interesting to study whether quantum effects of 
the breathing motion may give rise to fundamental changes in energy levels and phase 
structures of the skyrmion. These generic features including the phase structure 
seem to show close relations and striking similarities to the baryonic matter in the flat 
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space as known on classical level. 

In this paper we investigate the phase structure of a quantized chiral soliton on S^, 
where the effects of the quantization of breathing and rotational modes are taken into 
account on a basis of a family of trial functions for the profile function of the hegdehog 
ansatz.^)''^^ This approach is also a natural extension of the analysis in Ref. 6) for the 
skyrmion to those for skyrmion system. 

This paper is organized as follows. In Sec. II the standard skyrmion on is quan- 
tized by taking account of the breathing mode in addition to the spin-isospin rotation. 
In Sec. Ill on the basis of a family of trial functions for the profile functions of the 
skyrmion we obtain explicit expressions of the Lagrangian and the Schrodinger equation 
for the spin-isospin rotation and the breathing motion. In Sec. IV we investigate the 
effects of quantization of the breathing mode by solving numerically the Schrodinger 
equation for the / = J = state of the skyrmion. In Sec. V we study the phase 
structure of the quantized skyrmion obtained by solving the Schrodinger equation for 
spin-isospin states I — J — 1/2, 3/2 and 5/2. The final section is devoted to discussion 
and conclusions. 



§2. Breathing Motion of Skyrmion 

In the standard Skyrme model^^'^^ on S^, the Lagrangian of a chiral SU{2)l x 
SU {2)ji nonlinear cr model is given in the massless pion limit as 

La =\fl'Y^{9^''d^U{x)d.U\x)), (1) 

where U{x) = exp{2i7r{x)/ /t^) is the element of SU{2) group, 7r{x) = t^'^{x)t°'/2 de- 
scribes the pion field, /^r = 93MeV is the pion decay constant. The £sk term, called 



Skyrme term/) is introduced to stabilize the soliton solutions where e is a free param- 
eter. Here g^^ is the inverse of metric tensor g^y in a curved space where g^^ = (5o^. 
Using the standard spherical coordinate i = {u>, 9, (f)) on 5'^(p), the metric tensor gij on 
the 3-sphere of radius p are given as 

Qij = -Sijgi, guj = p"^, ge = sin^ a;, g^ = sin^ u sin^ 6. (2) 

The hedgehog ansatz for the classical static skyrmion-*^)'^) is given by 

Uq{x) = exp[zx • tF(lij)] 

X =ex sin 9 cos (f) + By sin 9sm(j) + cos 9, (3) 

where F{0) = tt and F(7r) = for the baryon number B = 1 sector. Substitution of Uq 
into U of (1) as 

Mo = -J C{Uo{x)) dV, 

dV = s/gi^x = p^ sin^ to dwdcj) sin 9 d9, g = det^r^ , (4) 

leads to the static energy functional;^) '^^ 

Mo[F]=27rp du sin 9d9[M A + Msk], 
Jo Jo 

Ma =^{f'^ sin^ u + 2sin^ F), (5) 
A<SK=^{^ + 2(F'sinF)2}, 

where F means F{u>) and F' = dF{u>)/du>. The classical profile function Fd minimizes 
the functional Mo[F]. 

The quantization of the rotational and breathing modes is performed by the fol- 
lowing replacement in the Lagrangian (l):^-* 

Uix.t) = ^(t)exp[zx • TFioj^xmm^, (6) 
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where A{t) e SU (2) and x(^) are the dynamical variables of the spin-isospin rotation 
and the breathing motion respectively. It is convenient to express the breathing motion 
as F{u!,x{t)) = F{iJ,) by use of a dimensionlcss variable ii ~ tan(a;/2)/x(t) which assures 
invariant boundaries of the range of variable fj, for any x(t) as = at a; = and fx = oo 
at a; = TT. If the boundary values of a given trial profile function F(yLt) are taken as 
F{0) = TT and F{oo) = 0, therefore, the breathing profile function F{uj,x{t)) satisfies 
automatically the boundary conditions for the baryon number B = 1 sector on S^, i.e. 
F{(jj = 0,x(^)) = TT and F{(jj = tt, x(^)) = 0. It is worth noting that for the limit of the 
3-sphere radius p — > oo, /x is just equivalent to the dimensionlcss variable r/R{t) of the 
flat spaces where R{t) is regarded as a dynamical variable for the breathing motion 
of the skyrmion.^)"'^^ For the simultanious limit p —>■ oo and a; ^ 0, a radial distance 
r between a point (a;, 9,(f)) on and the origin of R^ is approximated to 2ptan(a;/2) 
which is just the radial distance between the point (uj, $, </>) and the north pole (a; = 0) 
of S^. Thus, if we identify the dynamical variable xi^) with R{t)/{2p), where R{t) is 
a dynamical variable of the breathing motion in the flat space R^,^^ The dimensionlcss 
variable ji = tan(a;/2)/x(t) behaves as 

CO 

2ptan- ^ 
lim u = lim — — - — = ^, , . (7) 

Now the proflle function F{ii) = F{r/R{t)) describes the breathing motion of a localized 
skyrmion at the origin of R^. These are the reason why we use the proflle function 
with jj, = tan(u;/2)/x(t). 

By substituting (6) to (1) we obtain the Lagrangian for the spin-isospin rotation 
and the breathing motion of the chiral solitons given by 

L = a(x)x' + A(x)TV[iit] _ v{x). (8) 

Here a(x) and A(x) are the x(^) dependent inertia of the breathing motion and the 
moment of inertia of the spin-isospin rotation, respectively, and are given by 



a(x) =2n r du;p [fy sin^ u; + ^ sin^ F] (^) ^ 
^0 6 "X 



A(X) =^ du;p' sin^ [l + -^{F'^ + ^)] sin^ F, 



(9) 



where F = F(a;,x(^)) and F' — dF {u , x{t)) / du . It is noted that dF/dx is rewritten as 
— {li/x)dF{ii)/dii by use of an expression F{uj,x{t)) = -^(a*) with = tan(a;/2)/x(t). 
The V{x) = Mo[F] is the potential of the breathing motion obtained by replacement of 
F(a;)by F(a;,x(t))in (5). 

According to the standard procedure of the quantization based on the Lagrangian 
(8), after separation of the rotational mode, we obtain the Schrodinger equation for the 
wave function of the breathing mode, i/^j{x) for a / = J spin-isospin eigenstate:^)'^^ 

t-sb *l? + "'^'i' + ^^^^ ^ ^ - ^^i'^^'^' = 

where Ej is an energy eigenvalue and 

2X(x) 2aix)' 
and prime means differentiation with respect to x- 



§3. Trial Profile Functions and Quantized Skyrmion 



In order to obtain the quantized skyrmion solutions on S^, we have to solve the 
Schrodinger equation of the breathing motion (10) for spin-isospin eigen states varying 
the profile function which gives the most minimum eigenvalue Ej. To investigate 
for the wide range of variation of the profile function we adopt the variational method 
in terms of a family of trial functions F{jj,; C) with a parameter C {C > 0), which is 
expressed by an integral form as 

poo 7 
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where 

dx 

The trial profile functions C) satisfy the boundary condition F(0; C) = tt corre- 
sponding to a soliton of the baryon number 1 and the asymptotic condition for finite 
energy C) ^ for — > oo. C) with finite C has the asymptotic falloff oiXj y? 
for large [i which is required by the Euler-Lagrange equation for the profile function in 
the chiral limit regardless of the stabilizer. This asymptotic falloff also agrees with that 
of the static meson theory in the massless pion limit. 

In the usual skyrmion the function with C = 1, F{\i\C — 1), is just the 
same as the one derived by Atiyah and Manton^^ from the instanton solution and it is 
the nearest trial function to the classical numerical solution Fci{iJ,) among the family 
F{fi; C). It gives a static energy only 0.9% above that of numerical solution.^) The case 
with C = 2 is a trial function of arccosine form introduced by Igarashi, Otsu and two of 
the present authors/^^ which gives also a very good approximation of the exact static 
profile function Fci{iJ,) and gives a static energy only 1.0% above that of the numerical 
solution. 

By use of the trial profile function (11) we obtain the inertias ^(x) and A(x) from (9) 
and the potential V{x) = Mo[F] from (5). These inertias and potential are dependent 
of the parameter C. Varying this parameter C we solve the Schrodinger equations (10) 
for the breathing motion with definite spin and isospin I = J. 

In the limit of the radius p — > cxd the skyrmion tends to a usual skyrmion. 
In this case, when the Skyrme term contribution is small (i.e. when the Skyrme term 
constant e > 5.0 ), there appear the " frozen states" in which all the freedoms of the 
spin-isospin rotation and the breathing motion are completely frozen. This is because 
the state with infinite inertias of the breathing motion and with infinite moment of 
inertia of spin-isospin rotation, so that the contributions from the kinetic energies of 
both modes become zero and the skyrmion energy comes only from the potential 
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energy of the breating mode. The profile function in the frozen state is given by 

F{n; C = oo) which gives the most minimum eigenenergy E^o in the region e > 5.0. In 
this region of e the classical profile function Fci{n) or the trial function F{iJ,;C — 1) 
gives local maximum eigenenergies. 

On the other hand when the Skyrme term constant e is small the frozen state does 
not appear for the small spin-isospin states and the classical profile function which is 
approximated by F{ii; C = 1) gives most minimum eigenenergies. 

We confirmed numerically that in the case of skyrmion a similar situation occurs 
as in the skyrmion. The eigenenergies of skyrmion with the definite spin-isospin 
I = J become minimum at C = 1 or at C = oo depending on the parameter e. In the 
following, however, we calculate numerically for two special values of parameter C, i.e. 
C = 2 and C = oo. The value C = 2 is a substitute of C = 1. The reason for solving 
with C = 2 instead of C = 1 lies on the fact that we can obtain analytically integrated 
results of the potential V{x) and inertias a{x) and A(x) for C = 2 as well as C = oo 
but not for C — 1. The difference between the eigenenergies for C — 2 and C = 1 is 
expected within 1% as in the case of the skyrmion. 

[i) Breathing motion, inertias and potential for C = oo 

The tial profile function for C = oo 

2 

C = oo) = 7r(l arctan//), (12) 

is rewritten as ^ 

^^^2 F(ii;C = oo) 

which is known as a conformal map.^^'^-* 

Now we investigate the breathing motions of the profile function and baryon density 
given by the conformal map: 

F{u}, xit)) = 2 arccot /i. (12") 
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In Fig. 1(a) we show the breathing profile functions for various values of x(*)) 
X = 1/5, 1/3, 1, 3 and 5. The solid line represents the conformal map profile functions. 
The diagonal line represents the identity map Fiuj, 1) = u). It is noted that the peak of 
the profiles x(t)) at uj = shrinks or spreads out as the values x(^) decreasing or 
increasing from 1. 

The baryon density for the x(^) at the radial angle u on is defined as 



where is the baryon number density. The solid line in Fig. 1(b) represents the 

baryon number density given by the conformal map for various x(^) values. The dis- 
tributions of the baryon number density for the conformal map profile function have a 
exact symmetry. It is easy to prove by using the explicit form for F(a;,x(^)) of (12") 
that the effect of an inversion x ^ 1/x in (13) is equal to the one of a transformation 
a; — TT — u; by which the north semisphere of the is transformed on to the south 
semisphere. Thus, the peaked shapes of the baryon density around the north pole for 
X = 1/2 and 1/3 are exactly mirror symmetric to those of the corresponding one around 
the south pole for x = 2 and 3, respectively. The identity map Fioo, 1) = u gives com- 
pletely fiat baryon density over the whole 3-sphere as shown by the solid line of x = 1 in 
Fig. 1(b), i.e. 27r^j°(a;)p^ = 1. Thus, the breathing motion expressed by the dynamical 
variable x(^) means that distributions of the baryon density fiip-fiops between a higher 
peak at the north pole (e.g. x = 1/2) through completely fiat one (x = 1) and to higher 
one at the south pole (x = 2), and vice versa, as shown in Fig. 1(b). 

Here, by substituting (12) to (9), we obtain integrated expressions of the inertias 
a(x) and A(x) and the potential V{x) as follows. 



sin'' F{uj,x{t))dF{u,x{t)) 



(13) 



sin^ u) dijj 



a(x) 



47r2 



A(X) 

X^ 



A(X) 



IGtt^ 



A(X), 
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X+-+2 ^ 
X 

A(X) = 1 (3L^ + L{x+-+ 4)). 

(X + - + 2)2 >^ 
X 

where L — ef-j^p is the dimensionless radius and wiU be use instead of the radius p 
of the 3-sphere in the foUowing. 



We notice that under an inversion; 



X^K (15) 



there exists a symmetry as follows; 



A(x)=A(i), (16) 

Therefore for the conformal map profile function (12), the Lagrangian (8) and the 
Schrodinger equation (10) are invariant under the inversion (15). Hereafter, we call this 
invariance as the inversion symmetry. In Fig. 2(a), the potential V{x) of the conformal 
map case (12) is shown for varying the radius L. The potential exhibits inversion 
symmetry. In the region L < \pl the potential V^(x) has single valley symmetric for the 
variable Inx and the potential V^(x) takes the absolute minimum value at x = 1 for fixed 
value of L. This conformal map with x = 1 called identity map realizes a homogeneous 
distribution of the baryon number density over the whole 3-sphere as shown in Fig. 1(b). 
In these L region, the minimum of potential for the conformal map profile function is 
given by 

F„,i„ (conformal) = ^^(L + ^), L < v^, (17) 

e L 

It is noted that this really gives the energy of exact solution, i.e. the static classical 
energy Ed-, as shown in Ref. 3 and 4. Here, it is also worth noting that the theoretical 
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lower bound energy Gtt^/^/c is realized at x = 1 ^^d L = 1,^)'^^ i.e. the absolute 
minimum of the potential valley in Fig. 2(a). 

In the region L > v^, the potential V{x) behaves as the double well potential of the 
variable Inx for a fixed L and has local maxima at % = 1 and two local minima at x = 
Xim = V2L-1-^2{L^-V2L) < 1 and at X = 1/X«m = ^/2L-l + ^f2{L^^W^) > l- 
The minimum value of potential energy is given as a function of L by^-*'^-* 

yn,in(conformal) = ^^^(2a/2 - -), L> V2. (18) 

e L 

This minimum value is considerably larger than those of the exact solution^^'^\ which 
is approximated well by the solution for the arccosine profile function with C = 2 given 
in the following. 

(ii) Breathing motion, inertias and potential for C = 2 

The trial profile function for C = 2 in (11) is written by arccosine form; 

Fifi- C = 2) = arccos (l - ^J^J^)- (19) 

The dotted line in Fig. 1(a) show the breathing profile functions of the arccosine 
form for values of x(^) = 1/5, 1/3, 1, 3 and 5. The dotted lines of Fig. 1(b) display 
the corresponding baryon number densities. The distributions of the baryon number 
density has no symmetry under the transformation x 1/x a; — > tt — a; in contrast 
to those for the conformal map profile function. The baryon density has a peak at the 
north pole for x =1/3, 1/2 and 1, but has a bump and dip around the south pole for 
X = 3 and 2 (> 1.6) as shown by dashed line in Fig. 1(b). 

Now, we obtain analytically integrated expressions by substituting (19) to (9) for 
the inertias a(x) and A(x) and potential V{x) as follows: 

Kx)=:^{L'A,ix) + LB,ix)), (20) 

e L 
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where 



^i(x) 



8x(32 - 20v^+ (8 + ^pi)x + Gx") 



5i(x) 



V^(l + X)^(l + V2x)=^ 
3 + 18x + 42x2 + 42x^ + 7x^ 



B2U) 



2X(1 + X)^ 
X'(7 + I2x) 

X(157 + 1256x + 4251x^ + 7632x^ + 7243x^ + 3176x^ + 397x^) 



(21) 



128(1 + x) 



AM 



8x(3 - 2v^ + (2 - V2)x) ^ x(3 + 12x + llx') 



(1 + X)2(l + V2x) 2(l + x)4 

23(llx+-) 



53 (x) 



512 



In this arccosine profile function case there exists no special symmetry in contrast with 
the conformal map case. As shown in Fig. 2(b) , in the region L < \f2 the potential 
y(x) has also a single well but with no special symmetries and has a minimum at 
X ~ 1.6. The minimum of the potential, however, is considerably larger than that of 
the conformal solution given by (17) in this small L region. This means that in classical 
level, the stable soliton solution for L < -\/2 is given by the profile function of the 
identity map (the conformal map with x = !)• 

In the region L > a/2, the V{x) behaves as a double well potential which has a local 
maximum at x ~ 1-6 (saddle point) and two local minima with two difi^erent depths. 
The smaller one of these two minima exists in the region of x < 1-6. This minimum 
value given by the arccosine profile function (19) gives a very good approximation of the 
exact solution within about 1% deviation. Thus, for L > y/2 the exact classical solution 
given by the Fd is approximated well by the arccosine profile function F(//; C = 2). 
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§4. Effects of Quantization of Breathing mode 

In order to examine the effects of quantization of the breathing mode we study 
here the solutions of the Schrodinger equation for the breathing motion of the / = J = 
state. 

In Fig. 3 we display the arbitrarily normalized wave functions V'J=o(x) ^ region 
< X ^ 1 as a typical example of the solutions of conformal map with C — oo. The 
solid line and the dashed line represent the symmetric state and anti-symmetric state, 
respectively, for L = 1.3, at which the eigenenergy of the symmetric state takes the 
lowest value among all L region. The value is -E/=j=o = 870Mev for e = 7. The 
potentials V{x) ^or L — 1.3 and L = 10 with e — 7 also shown by dotted lines in 
Fig. 3. The potentials and the wave functions in the region x ^ 1 ^.re given from those 
in < X < 1 by F(l/x) and ±V'j=o(l/x); respectively. For L = 10 we get degenerate 
states as shown by dot-dashed line in Fig. 3. The absolute values of wave functions have 
sharp peaks at x = Xim ~ 0.03 and x = ~ 33 where the potential V{x) takes 

the minimum value as shown in Fig. 3 for < x ^ 1- The probability distributions to 
find the x is so localized at x — Xim and x — ^/xim, therefore, the quantum tunneling 
effect between these two regions is considered to be negligible. 

In Fig. 4, we show calculated results of eigenenergy £^/=j=o as a function of L for 
several values of the Skyrme term constant e. Solid and dotted lines represent those of 
the symmetric and anti-symmetric solutions, respectively, for C = oo (conformal map). 
Two lines of eigenenergies are completely degenerate for large L-region. Dot-dashed 
lines represent eigenvalues -E'j=7=o for C = 2 (arccosine form). For e = 7, we display 
also in Fig. 3 by dashed lines the minimum values of the potential yinin( conformal) 
given by (17) and (18) for the conformal map profile function with e = 7. Therefore, 
the difference between the solid line and the dashed line represents contribution from 
the kinetic energy of breathing motion. 

From these results for / = J = the effects of quantization of the breathing mode 
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are summarized as the foUowings: 

(i) In classical level, there exist two phases.^) One is the phase where both baryon 
number and energy densities are well localized at the north or south pole of in 
the region L > v^, which corresponds to the skyrmion with usual profile function 
numerically solved. The other is the phase where these densities are homogeneously 
distributed in the the region L < y/2. In this case the profile function is given by that 
of the identity map. In the quantum level, the quantum breathing motion in terms of the 
profile function for conformal map on means quantum fiip-fiop between a localized 
distribution at north pole (x < 1) and a localized one at the south pole (x > 1) through 
a uniform distribution (identity map with x = !)• 

(ii) As shown in the previous section and also in Fig. 2(a), the quantum breathing 
motion is described by a single- well potential in the region L < \/2, but it is described 
by a double-well potential in the region L > \pl. In the conformal map case the height 
of the potential barrier at x = 1 from the bottom of potential increases from to oo as L 
increases from -\/2 to oo. Thus, the ground state wave function changes smoothly from 
the one which has a peak at x = 1 to the one with double peaks at around potential 
minimum points xim and 1/xim as L increasing from to oo. 

(iii) The one of new features of quantization of the breathing mode is the appearance 
of pairs of eigenstates; symmetric and anti-symmetric under the inversion (15) for the 
conformal map profile function. In the region of small L less than ~ 2 the splitting 
between these two energy levels is rather large due to the lack of energy barrier between 
small and large x regions. On the other hand, in the region of large L larger than ~ 3 
this pair of states has an almost degenerate energy levels. This means the potential 
barrier around x = 1 become high enogh to forbid the penetration between the small 
and large regions of x- As for the arccosine case, the potential V{x) shows a single- 
well or double-well behaviors depending on the sphere size L but does not show any 
symmetry such as given by (15). Therefore, the eigenenergies for the arccosine case 
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have no degeneracy. In Fig. 4 only lower energy solutions are shown. 

(iv) For small Skyrme term constant e and in the large L region the arccosine profile 
function gives lower energy than the conformal map case. For e = 3, 4 and 5 the 
arccosine case gives lower eigenenergy solution in a region L > 3. For e = 7 the energy 
of the arccosine case is slightly lower than the one of the conformal map case in a region 
3 < L < 20, but the conformal map case gives the lowest energy solution for L > 20 . 
For the limit L ^ oo we get the same results as Ref. 6 on i?^, i.e. for values of e larger 
than a critical value eo = 5.6 all minimum energy values of spin-isospin J states, Ej, 
degenerate to a value limL — > oo£'c=oo(e) = 7.788/e GeV of the conformal map case. 

(v) As for the classical skyrmion, the lowest energy is always realized by the radius 
L = 1 for any e. For the quantized skyrmion, however, the radius L giving the lowest 
energy E'q-^'^ increases slightly as e increase. That is, the lowest energy is £^q**" = 
1.902GeV at L = 1.1 for e = 3, E^^"^ = 1.456GeV at L - 1.1 for e = 4, E^'"" = 
1.185GeV at L = 1.2 for e = 5 and E^'"" = 0.870GeV at L = 1.3 for e = 7, respectively. 



§5. Phase Structure of Skyrmion 

Now we discuss the quantized skyrmion obtained by solving the Schrodinger 
equation with I = J — 1/2, 3/2 and 5/2. In Fig. 5(a) ~ (d) we show the L-dependence 
of eigenvalues Ej obtained for J = 1/2, 3/2 and 5/2 and for the Skyrme term constant 
e =3, 4, 5 and 7. These eigenenergies of the symmetric and anti-symmetric states 
for conformal map profile function are represented respectively by the solid and dotted 
linesand almost degenerate in the large L region. The dot-dashed lines in Fig. 5 represent 
the lowest eigenenergies for the arccosine profile function. 

For e = 3, as shown in Fig. 5(a), the arccosine profile function realizes the lowest 
energy solutions for all spin-isospin states I = J =1/2, 3/2 and 5/2 in the region L > 2. 
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In the region L < 2 the lowest energies are obtained by the conformal map profile 
function for all spin-isospin states. Then a phase transition occurs at L ~ 2. 

Fig. 5(b) exhibits the L-dependence of eigeneneries for e =4. In this case, for spin- 
isospin states 7 = J = 1/2 and 3/2, there occurs a phase transition at L ~ 2.5 i.e. the 
conformal and the arccosine profile functions give lower energies in the region L < 2.5 
and L > 2.5, respectively. On the other hand, for / = J = 5/2 the conformal map 
profile function always gives the lowest energy. 

For e = 5, as shown in Fig. 5(c), in a region L > 2.5, the arccosine profile function 
gives the lowest energy solution for / = J =1/2 but the conformal map profile function 
gives the lowest one for / = 1/2 in the region L < 2.5 and for / = J = 3/2 and 5/2 in 
all region of L > 0. Then the phase transition occurs at L ~ 2.5 only for J = 1/2, but 
not for J = 3/2 and 5/2. 

Fig. 5(d) shows for the case e = 7. For all spin-isospin states with J = 1/2,3/2 
and 5/2, the conformal map profile function always gives the lowest energy solutions in 
whole L region. 

We have considered hitherto the case where the radius of is fixed a priori by 
some condition given from outside. Now we turn to a rather academic or fictitious 
problem where the isolated skyrmion is putting on the vacuum and its radius p or L 
is determined by itself. In this case we have to regard the sphere radius L as a dynamical 
variable and to find minimum eigenenergies by solving Shrodinger equation for both L 
and X- If W6 assume that the radius L is an adiabatically variable of the skyrmion, 
then we have to search the lowest eigenenergies varying the radius L. In Table I we list 
the local minimum energies Ej in the small L region of symmetric states and the values 
of L which give the minimum energies for e = 3, 4, 5 and 7. Fig. 6 shows e-dependence 
of these local minimum eigenenergies Ej for J = 0, 1/2, 3/2 and 5/2. As shown by solid 
lines of Fig. 5(a) ~ 5(d), in the limit of L = oo the eigenenergies of the any spin-isospin 
states with conformal profile function converge to the same energy E^o = 7.788/e GeV 
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which is the energy of the frozen state of the skyrmion. The classical skyrmion 
energy E^i = QT^'^fj^/e and E^o are also shown by dashed and dot-dashed lines. If the 
5'"^ skyrmion with definite spin-isospin which is represented by the solid lines in Fig. 6 
above the dot-dashed line of £^00 the radius of skyrmion increase infinitly and the 
skyrmion given by the conformal profile functions disperses. The energy of / = J = 1/2 
state is always smaller than £^00 in the region e. For e > 5 the S"^ skyrmion with 
I — J — disperses while for e < 5 the / = J = 3/2 skyrmion given by the conformal 
map profile function has most minimum energy. The I = J = 5/2 skyrmion of 
conformal map becomes most minimum energy state for e < 4.1. Occurrence of infinite 
dispersion of the quantized skyrmion is recognized by the e and L dependence of the 
kinetic energies of the rotational and breathing modes and the potential of breathing 
mode. As given in (14) and (20) both inertias a(x) and A(x) are proportional to 
and increase as L increases. Then the kinetic energies of the rotational and breathing 
modes to and decrease as the L increases. On the other hand the potential energy 
of the breathing mode is proportional to e~^. Therefore in the large e region where the 
contributions of kinetic energies dominate, increase of L reduces total eigenenergy and 
eventually the skyrmion disperses as L — > 00. 

As listed in Table I, the values of minimum points L which give the lowest eigenen- 
ergies Ej are shifted toward larger values from that of the classical one, L = 1. This 
is also explained by the quantum effects of the kinetic energies of spin-isospin rotation 
and breathing motion. 



§6. Discussion and Conclusions 

In this paper we have investigated the quantum effects associated with rota- 
tion and breathing motions of a skyrmion making use of a family of trial profile 
function F(//; C) with fj, = tan(a;/2)/x(t) where x(^) represents the dynamical variable 
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of breathing motion of skyrmion. 

In the classical level, for L < we have a skyrmion whose profile func- 
tion is given by identity map F{u!) = u> and its energy density is distributed homoge- 
neously inside S^. In this case the baryon number density is l/2n'^p^ = e'^/^/27r^L^ > 
e^/^/4v^7r^. If the Skyrme term constant e is identified to the coupling constant 
5'p7r7r=5.8, then the baryon number density is larger than 0.37fm~^ which corresponds 
to high baryon number density matter and considered as deconfined state. ^^'^^ As shown 
in the preceding sections, in the region L < \f2 where the potential V{'x) of the breath- 
ing motion has single valley at x = 1, the quantiezed skyrmion is realized by the 
conformal map profile function. The ground state wave function of breathing motion 
of the / = J = /S^ skyrmion is shown by a full line in Fig. 3 which is symmetric 
under inverse transformation (15). As in the classical level this state is considered as 
the deconfined phase of baryonic matter, then the breathing motion obtained here can 
be regarded as a part of the breathing motion of the baryonic matter. 

In the region L > \/2, for most region of parameters e and L and for the states 
with lower spin-isospin, we obtain a localized skyrmion in the classical level whose profile 
function is given by arccosine form (19) and considered as the confined state. ^^'^^ The 
arccosine profile function is an approximation of the Atiyah-Manton profile function 
given by (11) with C = 1. In this case the potential V{x) has two valleys as shown 
in Fig. 2(b) and the one at x smaller than 1.6 is deeper than the one at x larger than 
1.6. Then the wave function of the lowest energy is localized at very small x and the 
quantized skyrmion is more sharply confined at the north pole of S^. In the quantum 
level, however, there occurs for a certain region of the parameters e and L and for the 
spin-isospin states that the conformal profile function gives the lowest energy state of 
skyrmion. We take L = 10 as an example. The I = J = 1/2 skyrmion is given 
by the arccosine profile function in the region e < 6.5. Only in the region e > 6.5 the 
conformal profile function gives stable skyrmion. The skyrmion with / = J = 3/2 
is realized by the arccosine and the conformal profile functions in the regions e < 4.8 and 
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e > 4.9, respectively. For I = J = 5/2 skyrmion state the conformal profile function 
gives the lowest energy in the region e > 4 and the arccosine profile function does in the 
region e < 4. In general in the region of small e corresponding to strong stabilizer, the 
lower spin-isospin states are given by the arccosine (or Atiyah-Manton) form of profile 
function and localized at the north pole of for large L. When the Skyrme term 
constant e increases, i.e. the contribution of stabilizer becomes weak, the skyrmions 
change to the ones given by the conformal map corresponding to deconfined system 
from higher spin-isospin states and eventually all spin-isospin states. 

The conformal map profile function corresponds to the " frozen states" on the 
skyrmion where the inertias a{x) and A(x) of the breathing and rotational modes are 
infinite. However, since in the case of the skyrmion the value L is finite then the 
inertias a{x) and A(x) are also finite contraly in the case of skyrmion, the " frozen 
states" do not appear any more. 

In conclusion, there appears two types of breathing motions, one is the breathing 
motion of the profile function given by conformal map, and the other is those approx- 
imated by the arccosine profile function. For strong stabilizer the lowest eigenenergy 
state in large L region is given by the breathing motion of arccosine profile function. 
For the weak stabilizer the eigen state of (10) for the conformal map profile function 
always gives the lowest energy in all radius L, which continuously connects to the 
"frozen states" on R^ in the limit L — > oo. Howeve, in the small L region, the 
lowest eigen state is always given by the breathing motion of the conformal map profile 
function. 

Since the skyrmion approaches provide very good modeling of the periodic arrays 
of the cubic lattice skyrmion in fiat space R^ as shown by various studies in the classical 
level, above results are not specific ones of the skyrmion, but also essential features 
of the quantized lattice skyrmions in flat space R^. For example, the quantum breathing 
motions on around the uniform distribution (identity map with x = 1), correspond 
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to the quantum fluctuations of the baryon density around the homogeneous distribution 
with the half-skyrmion symmetry^) of the periodic array of the cubic lattice skyrmion 
in the flat space R^. Thus, it may be possible to flnd these remarkable features in the 
dense nuclear matter. 
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Figure Captions 

FIG. 1 (a) The breathing profile functions for x =1/5, 1/3, 1, 3 and 5. The sohd fine 
represents the conformal map profile functions and the dashed line represents the 
arccosine profile function, (b) The baryon density as a function of the radial angle 
uj for X =1/3, 1/2. 1, 2 and 3. The solid line (dashed line) represents the baryon 
density for the conformal map profile function (the arccosine profile function). 

FIG. 2 The potential V{xiL) as a function of two variable x ^-^d L; (a) the potential 
for the conformal map profile functions, (b) the potential for the arccosine profile 
function. 

FIG. 3 The potentials V{x) in GeV and the arbitrarily normalized wave functions "00 (x) of 
the spin-isospin state with I = J = for the Skyrme term constant e — 7 plotted 
against the variable x- The dotted line shows the x dependence of the potential 
V{x) in GeV for L = 1.3 or for L = 10. The solid line (dashed line) represents 
symmetric (anti-symmetric) wavefunctions V'o(x) for L = 1.3, and dot-dashed line 
represents a wavefunction for the almost degenerate state for L — 10. 

FIG. 4 The eigenenergy £^/=j=o in GeV of three types of eigen states as a functions of L for 
various Skyrme term constants e; the symmetric states (solid line), anti-symmetric 
states (doted line) and almost degenerate states (solid line) for conformal map 
profile function, and almost degenerate states (dot-dashed line) for the arccosine 
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profile functions. The dashed fine represents the minimum value of the potential 
as a function of L for the conformal map profile function with e = 7. 

FIG. 5 The eigenenergy Ej^j in GeV of the spin-isospin states with / = J = 1/2, 3/2 and 
5/2 as a function of L for various Skyrme term constants e; (a) e = 3, (b) e = 4, 
(c) e = 5 and (d) e — 7. For the conformal map profile functions, the energy of the 
symmetric state , anti-symmetric state and almost degenerate state are represented 
by the solid line, the dotted line and the solid line, respectively. The dot-dashed 
line represents the lowest eigenenergy for the arccosine profile functions. 

FIG. 6 The e-dependence of the lowest eigenenergies Ej for J = 0, 1/2, 3/2 and 5/2. The 
classical skyrmion energy E^ = Qir'^fT^/e and E^o are also shown by dotted and 
dot-dashed lines. 
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TABLE I. The lowest energies EJ'"" in MeV. 



e 


L 


pmin 


L 


^1/2 


L 


-^3/2 


L 


TTTnin 

^5/2 






3 


1.1 


1902 


1.1 


1931 


1.2 


2035 


1.4 


2180 


1836 


2596 


4 


1.1 


1456 


1.2 


1518 


1.4 


1712 


1.8 


1928 


1377 


1947 


5 


1.2 


1185 


1.4 


1294 


1.8 


1566t 


2.5 


1837t 


1101 


1558 


7 


1.3 


870 


1.8 


1090 


2.7 


1455t 


3.5 


1750t 


787 


1113 



The state a f is attached to is not stable and will infinitely disperse when L is 
varied. 
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